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ABSTRACT 

We present spectroscopic confirmation of ten highly luminous (L ^ 2L*) Lyman alpha 
emitters in the redshift range 6.01 < z < 6.49 (nine galaxies and one AGN), initially 
drawn from a sample of fourteen z p hot ^ 6 Lyman break galaxies (LBGs) selected from 
an area of 0.25 square degrees within the UKIDSS Ultra-deep Survey (UDS). Overall, 
our high rate of spectroscopic confirmation 71%) and low rate of contamination 
provides a strong vindication of the photometric redshift analysis used to define the 
original sample. By considering star-formation rate estimates based on the Lya and 
UV continuum luminosity we conclude that our sample is consistent with a Lya escape 
fraction of ~ 25%. Moreover, after careful consideration of the potential uncertainties 
and biases, we find that 40% — 50% of our sample of L ^ 2L* galaxies at 6.0 < z < 6.5 
display strong Lya emission (rest-frame equivalent width ^ 25A), a fraction which is 
a factor of ~ 2 higher than previously reported for L ^ L* galaxies at z ~ 6. Our 
results suggest that, as the epoch of reionization is approached, it is plausible that 
the Lya emitter fraction amongst luminous (L ^ 2L*) LBGs shows a similarly sharp 
increase to that observed in their lower-luminosity (L L*) counterparts. 
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1 INTRODUCTION 

Improving our understanding of the earliest epochs of galaxy 
formation and evolution relies fundamentally on the ability 
to select clean samples of high-redshift galaxies, free from 
significant low-redshift contamination. Traditionally, sam- 
ples of high-redshift galaxies have been selected using one of 
two complementary photometric techniques. Firstly, Lyman- 
alpha emitters (LAEs) are selected from deep imaging using 
narrow-band filters centred on the redshifted Lyman-alpha 
emission line (e.g. Hu, Cowie Sz McMahon||1999 1. Alterna- 
tively, Lyman-break galaxies (LBGs) can be selected from 
deep broad-band photometry using the Lyman-break, or 
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"dropout", technique pioneered by |Steidel, Pettini fc Hamil-| 
~ 1 19951. 



ton 



Studies of high-redshift galaxies selected using both 
techniques have made rapid progress recently, thanks to 
large-area, red-sensitive, detectors on ground-based tele- 
scopes and the ultra-deep near-infrared imaging now pos- 
sible with WFC3/IR on-board the Hubble Space Telescope 
(HST). As a result, it is now possible to obtain large, statisti- 
cal, samples of LBGs/LAEs in the redshift interval 6 < z < 7 



from the ground (e.g. Yoshida et al. 2006 Ouchi et al 



|2008[ |2010| |McLure et al.||2009| and upwards of one hun- 
dred LBGs have now been identified in the redshift interval 



6.5 < z < 8.5 with WFC3/IR (e.g. Bouwens et al. 2010 



McLure et aLl|2010| |2011| |Finkelstein et al.||2010j )Tlndeed, 



recent results have demonstrated the power of these tech- 
niques with the spectroscopic confirmation of LBG-selected 
galaxies within the reionization epoch at 7.00 < z < 7.20 
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(Vanzella et al. 1120 11 Pentericci et al. 112011 Schenker et al 



2011 Ono et al. 20111. Although LAEs are simply a sub- 



set ol the LBG population, due to the different selection 
techniques employed, the study of these two high-redshift 
galaxy populations has proceeded largely independently. As 
a consequence, one of the key questions in the study of high- 
redshift galaxies is determining how LAEs and LBGs are re- 
lated and how they are connected to the galaxy populations 
identified at lower redshift. 

Over the last fifteen years substantial observational ef- 
fort has been invested in studying the population statistics 
of both LAEs and LBGs. As a result significant progress 
has been made towards an understanding of the luminos- 
ity function and clustering properties of the LAE and LBG 



strong Lya emitters between 6 < z < 7, though this is due 
mainly to changes in dust distributions, not a change in the 
neutral fraction of the IGM. However, they also conclude 
that all LBGs at z ~ 6 should display Lya in emissio n with 
EW>20A, in disagreement with the observations of Stark 



et al. (20111. Forero-Romero et al. (20111 are able to re- 



produce very well the observed fraction of EW>55A Lya 



emitters in the 4.5 < z < 6.0 sample from Stark et al. (2010), 



by requiring the escape fraction of Lya photons to decrease 
with increasing UV-luminosity. In addition, they are are also 
able to reproduce the observed drop in EW>25A Lya emit- 



ters above z > 6.3 observed by Schenker et al. (20111, by 



populations from z = 3 to z = 7 (e.g. Bouwens et al.||2007 



20T0l|2011||McLure et al.|2009|[20Tol|Ouchi et al.|2008||2010 



fitting to the observed evolution in the LAE luminosity func- 
tion at these redshifts. Neither of these studies allow for 
peculiar velocities (outflows/inflows) within the ISM which 
Dijkstra et al. ( 2011[) find to produce a large effect on the de- 



Reddy et al. 20081. While correlation length estimates for 
LAEs/LBGs are largely consistent (ro — 5 Mpc) and show 
little evidence for evolution, the LAE/LBG luminosity func- 
tions appear to show substantial differential evolution. Most 
studies now agree that the LBG luminosity function evolves 
strongly with redshift, with L* dimming by a factor of ^ 2 
between z = 4 and z = 6 (e.g. Bouwens et al.|2007 McLure 



et al.||2009| [201(5} |Yoshida et al.||2006[ ). In contrast, the ob 
served LAE luminosity function does not appear to evolve 
between z — 3 and z = 6 (e.g. Shimasa ku et al.|2 006 Ouchi 
et al.|2008 1, although the latest results indicate a moderate 
~ 30% drop in L* Lya between z = 5.7 and z = 6.6 (Ouchi 
et al.|2010| >. 



Studying the Lya emission properties of LBG-selected 
samples is a powerful method for improving our understand- 



tectability of Lya, with outflows from the galaxies boosting 
the observed EWs. 

Despite the large amounts of effort invested in spec- 
troscopically observing the properties of LBGs over recent 
years, one area of parameter space which has been relatively 
unexplored is that occupied by luminous LBGs at z ^ 6. The 
relatively poor statistics in this redshift-luminosity regime 
are simply a reflection of the relative rarity of L L* galax- 
ies at z ^ 6, combined with the small areas typically covered 
by survey fields with the deep, multi-wavelength, imaging 
necessary to reliably select such objects. In this paper we 
address this issue by presenting deep, red sensitive, optical 
spectroscopy of luminous (L 2L*) LBGs in the redshift 
interval 6.0 < z < 6.5, photometrically selected from an 
area of 0.25 sq. degrees within the UKIDSS (Lawrence et al. 



ing of this confusing picture (eg. Stark et al. 2010 Jiang 20071 Ultra-Deep Survey (UDS). Armed with deep spectro- 



et al. 20111. Although spectroscopic follow-up of the LBG 



population reveals Lya emission in objects covering the full 
span of UV luminosity, the strongest Lya emission (equiva- 
lent width > 100A) is only observed in fainter objects. This 
was first indicated by Shapley et al. (2003), who observed 



scopic observations, the combination of large area and deep 
optical/near-infrared imaging available in the UDS allows us 
to investigate three important issues. Firstly, we are able to 
investigate whether or not the photometrically selected LBG 
samples used to constrain the evolution of the bright end 



the mean rest-frame Lya equivalent width (EW) increasing 
at fainter magnitudes for z ~ 3 LBGs, and has since been 



noted by many authors over a wide redshift range (eg. Ando 



et al. 2006 Vanzella et al. 2009 Stark et al. 20101. Fur 



of the UV-selected galaxy luminosity function (e.g. McLure 
et al.||2009| ) are significantly contaminated by low-redshift 



thermore, studies of how the fraction of strong Lya emitters 
evolves as a function of redshift offer the prospect of pro- 
viding crucial information on the dust content of LBGs and, 
potentially, on the neutral fraction of the IGM. Recent work 
by |Stark et ah] ( |2010| |2011[ ) has shown that, at fixed UV 
luminosity, the fraction of LBGs showing Lya emission with 
EW > 25A increases by 55% between z — 4 and z — 6. More 
recently this analysis has been extended to spectroscopically 
targeted z ~ 7 samples ( Fontana et al.| 2010 Schenker et al 



2011||Pentericci et al.|2011||Ono et al.|2011[ ) with the results 
from different groups consistently showing a drop in the Lya 
fraction at z > 7, tentatively attributed to a change in the 
neutral fraction in the IGM. 

Recent theoretical studies find differing results when 
trying to explain the observed fractions of Lya emitting 



galaxies amongst LBG samples. Dayal & Ferrara (2011 \ are 



able to broadly reproduce the observed trend of increas- 
ing Lya emission with decreasing UV luminosity for an 
EW>55A cut at z ~ 6 from their cosmological SPH simula- 
tion of reionization. In agreement with recent observations, 



interlopers and/or active galactic nuclei (AGN). Secondly, 
we are able to place the best available constraints on the 
fraction of strong Lya emitters amongst luminous LBGs at 
z ^ 6. Finally, we are able to investigate whether previous 
reports of a very low fraction of strong Lya emitters amongst 
samples of luminous z ~ 6 LBGs have been biased due to 
contamination by low-redshift interlopers or small number 
statistics. 

The structure of the paper is as follows. In Section 2 we 
describe the initial sample selection and spectroscopic obser- 
vations before proceeding to present the spectra themselves. 
In Section 3 we describe our technique for accurately mea- 
suring the Lya equivalent widths (EWs) and investigate the 
star-formation rates (SFRs) and Lya photon escape frac- 
tion of our LBG sample, before estimating the fraction of 
luminous z — 6 LBGs which are strong Lya emitters. In 
Section 4 we compare our results to those from the recent 
literature and discuss the potential uncertainties and biases 
associated with estimating the Lya emitter fraction at high 
redshift. In Section 5 we present our conclusions. Through- 
out the paper we assume a cosmology with Ho = 70 km 
s -1 Mpc - \ £lm = 0.3, fi A = 0.7. All magnitudes are quoted 



Dayal & Ferrara ( 2011[ ) also find a decrease in the fraction of in the AB system ( Oke fc Gunn|[l983 \ 
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Figure 1. Spectra of the confirmed z > 6 objects showing the one-dimensional (ID) spectra to the left and the corresponding 
two-dimensional (2D) spectra on the right. All ID spectra are plotted on the same flux scale, except for the spectrum of the AGN 
(UUDS-J021627.8) which has been scaled by a factor of 0.5. 



2 IMAGING DATA, SAMPLE SELECTION 
AND SPECTROSCOPY 

The sample of z ^ 6 LBG candidates targeted for spectro- 
scopic follow-up is a sub-set of that originally selected by 
McLure et al. (2009), using photometric redshift analysis, 
to investigate the bright-end of the z = 6 galaxy luminosity 
function. Although a full description of the original selection 
process is provided in McLure et al. (2009), in this section 
we briefly review the most relevant details before proceeding 
to describe the spectroscopic observations. 



2.1 Imaging data and LBG selection 

The parent sample of z 6 LBG candidates was originally 
selected via a photometric redshift analysis, which exploited 
the deep optical and near-infrared imaging available within 
the area covered by the UKIDSS Ultra-Deep Survey (UDS). 
The UDS is the deepest of the five near-infrared surveys 
being undertaken at the UK infra-red telescope (UKIRT) 
which together comprise the UK infra-red Deep Sky Surveys 
(UKIDSS; Lawrence et al. 2007). 

The latest ESO-public data-release for the UDS (DR8) 
features near-infrared imaging over an area of 0.8 square 
degrees to 5a-depths of J = 24.9, H = 24.2 & K = 24.7 
(2"diameter apertures). Within the UDS field, complimen- 
tary deep Subaru optical imaging is available from the 
Subaru/XMM Deep Survey (SXDS), which provides opti- 
cal imaging to 5<r depths (2"diameter apertures) of B — 
27.9, V = 27.3, R = 27.2, z' = 27.2 & z' = 26.1 (Furusawa et 
al. 2008). 



The useful overlap region between the UDS near- 
infrared and SXDS optical imaging covers an area of ~ 0.65 
square degrees, and it was from this overlap region that 
McLure et al. (2009) used a photometric redshift analysis 
to select luminous (z' < 26) LBG candidates in the redshift 
interval 4.5 < z < 6.5. The original selection utilised the full 
information available from the SXDS optical bands and the 
UKIDSS DR1 J— and K— band near-infrared imaging. Our 
updated photometric redshift analysis, using the UKIDSS 
DR8 J, H and K-band imaging, is described in greater de- 
tail in Section 2.2.3. The parent sample for spectroscopic 
follow-up was comprised of those candidates with a high 
probability of lying at z ^ 6 based on the redshift proba- 
bility density function returned by the photometric redshift 
analysis. 



2.2 Spectroscopy 

The spectroscopic data analysed in this study were obtained 
between October 2007 and January 2011 with the FORS2 
spectrograph on the VLT as part of the systematic spec- 
troscopic follow-up of the UDS obtained through the ESO 
large programme ESO 180.A-0776 (UDSz; RI. O. Almaini). 
Given that full details of UDSz will be presented in Almaini 
et al. (2012, in preparation), only the most relevant details 
are provided here. 

The UDSz programme was allocated a total of 235 hours 
of observations, with 93 hours allocated for observations 
with the VIMOS spectrograph and 142 hours allocated for 
FORS2 observations. The primary science driver for UDSz 
was to obtain spectroscopic observations of a representative 
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Table 1. The spectroscopic sample of 2 > 6 UDS LBGs. Column 1 lists the UDS object IDs (incorporating the J2000 
coordinates). Column 2 gives the spectroscopic redshift and errors derived from the FORS2 spectra, taking the spectroscopic 
rcdshift from the peak of the Lya line. Columns 3 & 4 give the quality of the spectroscopic redshift (A-C, see text for 
details) and the photometric redshift derived in Section 2.3 respectively. The remaining columns list the z'-band and 2:921- 
NB observed magnitudes (corrected to total) plus errors, absolute UV magnitude derived from the 2921-NB photometry 
and the measured Lya fluxes and EWs (derived from spectra, see text for more details). 
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a Lya flux calibration tied to the z— band photometry using an average correction factor, see text for details. 

^ EW values potentially affected by Lya contributing to the NB photometry used to derive the UV continuum value. 



sample of if —band selected galaxies (K < 23) photometri- 
cally pre-selected to lie at redshift z v hot =S 1-0. These spec- 
troscopic observations were designed to support the primary 
science driver of the UDS survey, which is to study the as- 
sembly and evolution of massive galaxies at z ^ 1. Within 
this context, the large allocation of FORS2 time (20 MXU 
masks, 5 hours of integration each) was designed to pro- 
vide sensitive, red optical, spectra of massive galaxies in the 
redshift interval 1.0 < z < 1.5. 

Each of the 20 FORS2 masks covers an area of 6.8' x 
6.8' and typically allowed for the allocation of 30-35 science 
slits (based on a slit length of 8"). During the mask design 
process targets from the core sample of K < 23 galaxies 
were allocated with the highest priority. However, in order 
to fill each mask, several other samples of additional sci- 
ence targets were also allocated slits on a best-effort basis. 
The list of high-redshift LBG candidates from McLure et al. 
(2009) with a high probability of lying at z ^ 6 was included 
as one of the samples of additional science targets. A total 
of fourteen z ^ 6 LBG candidates were observed as part of 
UDSz, each receiving a total of 5-hours of on-source integra- 
tion with the GRS-300I grism (6000A< A < 10000A with 
R — 660). Data reduction was performed using a modified 
version of the FORS2 pipeline, details of which are provided 
in Pearce et al. (2011, in preparation). 

2.2.1 Spectra of confirmed 2^6 sources 

Of the fourteen galaxies for which spectra were obtained 
we are able to assign secure redshifts to eleven sources ac- 
cording to Lya line detections, placing them at z > 6. The 
ID spectra for the spectroscopically confirmed z > 6 Lya 
emitters are displayed in Figure [T] and the measured line 
fluxes and EWs are reported in Table [I] Many of the spec- 
tra show continuum red-wards of the emission line as well as 
the distinctive asymmetrical line profile indicative of Lya. 
Each redshift is assigned a quality flag, defined as in |Vanzella| 
et al. ( 2009[ ), with values of A (unambiguous detection), B 
(likely detection) and C (uncertain detection). It should be 
noted that the redshifts for the 11 confirmed 2 > 6 sources 



are measured from the position of peak flux within the Lya 
line, and are hence susceptible to shifts in the peak position 
due any velocity structure of the neutral hydrogen within 
the galaxy or the surrounding IGM. The average velocity 
shift measured for Lya in lower redsh i ft samples is of orde r 
+350 km s _1 |Adelberger et al.||2003| |Shapley et aT|2003| ), 
since the Lya photons are more likely to escape the galaxy 
when they are scattered off neutral Hydrogen that is travel- 
ling away from us, giving the photons a Doppler shift that 
takes them off resonance. 

One of our targets (UUDSJ021627.80-045534.2) is a 
low-luminosity active galactic nucleus (AGN) at a redshift 
of z = 6.01 (Lya FWHM ~ 1600 km s" 1 ). This object was 
initially identified as a potential ultra-luminous z ~ 6 LBG 
candidate by McLure et al. (2006) and was subsequently 
targeted as a high-redshift quasar candidate by Willott et 
al. (2009). In Willott et al. (2009) this source is referred to 
as CFHQS J021627-045534, and was confirmed as a low- 
luminosity quasar at z — 6.01 based on a four-hour spectrum 
obtained with the GMOS spectrograph on Gemini. However, 
it is important to note that CFHQS J021627-045534 was 
the only z ~ 6 AGN identified by Willott et al. (2010) over a 
search area of 4.47 sq. degrees, down to a limit of z' = 24.5. 
Based on this fact, and the steepness of the z ~ 6 UV- 
selected galaxy luminosity function, it seems reasonable to 
conclude that the luminous, type I AGN contamination of 
26 ^ z' 25 LBG samples at z ~ 6 is negligible. 

Based on the optical spectra alone, we cannot entirely 
rule out the possibility that the observed Lya emission could 
have some contribution from a narrow-line AGN component. 
However, we note that a recent study of UV-selected AGN 
at z ~ 2 — 3 by |Hainline et alT ( 2011 1 indicates that narrow- 
line AGN are on average much redder than non-AGN LBGs 
(/? = -0.314 compared to (3 = -1.49) and show the NV 
line at \ r est = 1240A in emission. Given that our sample 
have broad-band SEDs which are not well reproduced by 
a standard Type-2 AGN template (see Figure [SJ, are not 
detected at 24/xm and do not display NV emission lines in 
their FORS2 spectra, we conclude that the UV continuum 
in these objects is dominated by starlight. 
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Figure 2. SED fits to the spectroscopically confirmed z ^ 6 LBGs showing the BV RizJ H K photometry (red data points with error 
bars), the best fitting template with a fixed spectroscopic redshift and including the measured Lya contribution (black line and blue 
data points) and the best-fitting template assuming no Lycr flux and leaving redshift as a free parameter (grey lines). The best-fitting 
SED templates at low redshift (z p h t ^ 2) are plotted as dashed grey lines. The inset panels show the distribution of x 2 versus redshift. 



2.2.2 Spectra of unconfirmed sources 

Four of the z ^ 6 LBG candidates which were spectro- 
scopically observed did not produce spectra which allow 
us to assign a robust redshift. One of these four objects 
(UUDS_J021816. 33-051116. 6) does display a probable Lyo? 
emission line but is assigned a quality flag of C since the 
emission line sits within a region of poorly subtracted sky- 
line residuals. One further object does not display a promi- 
nent Lya emission line, but does display a continuum break 
consistent with a redshift of z = 5.83. However, the signal- 
to-noise in the final FORS2 spectrum is too low to be con- 
fident of the redshift. 

The spectra of the remaining two candidates show only 
a very faint continuum, with a signal-to- noise ratio too low 
to reliably distinguish between a Lyman break and low- 
redshift interloper. We note that three of the four uncon- 
firmed candidates lie on the same FORS2 mask. However, 
although this mask is slightly poorer than average, in terms 
of redshift completeness for primary targets at 1 < z < 1.5, 



there is no particular reason to suspect that the spectra of 
these three high-redshift candidates have been unduly af- 
fected. 



2.3 Photometric redshift selection and 
spectroscopic completeness 

From the sample of fourteen z ^ 6 LBG candidates which 
were selected for spectroscopic follow-up, a total of ten 
objects have been robustly spectroscopically confirmed as 
z ^ 6 objects. As a consequence, regarding the other four 
spectra as contaminants, the lower limit to the spectroscopic 
completeness for our sample is 71%. However, as described 
above, two of the remaining spectra display strong evidence 
of being genuine high-redshift galaxies at z 5.8 and have 
only failed to make our robust sample due to insufficient 
signal-to-noise. If we regard these two objects as also be- 
ing confirmed z ~ 6 objects, then the spectroscopic com- 
pleteness of our sample rises to 86%. Irrespective of the na- 
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Figure 3. SED fits for the three unconfirmed objects and UUDS-J021816. 33-051116. 6, which is assigned a quality 'C emission line 
redshift (see text) and is not treated as a secure Lya detection. For the three objects without an assigned spectroscopic redshift, only 
the primary and secondary photometric redshift solutions are shown as grey solid and dashed lines respectively. For these objects the 
model fluxes are calculated from the primary redshift solution (blue points). 



ture of the two uncertain candidates, the high spectroscopic 
completeness provides a strong vindication of the original 
method of selecting high-redshift targets on the basis of an 
SED-fitting photometric redshift analysis. 

To illustrate this point, in Fig. 2 we show the latest 
photometry available for each of the spectroscopically tar- 
geted z ^ 6 objects, along with the results of our updated 
photometric redshift analysis. For each object the observed 
photometry is plotted along with the best-fitting SED model 
returned when redshift is kept as a free parameter, as well 
as the best-fitting model using the fixed spectroscopic red- 
shift and adding a Lya emission line with the correct EW. 
The photometric redshifts were calculated using LePhare 
photometric redshift code ( Ilbert et al.||2006 I. The galaxy 



templates were made using the Bruzual & Chariot ( 2003 1 



stellar population synthesis models, with exponentially de- 
creasing star formation rates with e-folding times in the 
range 0.1 Gyr < r < 30 Gyr and a Chabrier (20031 ini- 



tial mass function (IMF). Reddening by dust is added using 
the |Calzetti et al.| ( [2000[ ) extinction law with E(B-V) values 
ranging from 0.0 to 0.5. For the fitting with fixed redshift, 
Lya is added to the models after extinction is applied, using 
the flux from the template convolved with the narrow-band 
filter as the estimate of the continuum to translate the EW 
into the Lya flux. 

In Fig. 3 the insets showing the distribution of \ 2 versus 
photometric redshift clearly demonstrate, as expected, that 
the primary photometric redshift solutions are all at z>6. 
Crucially, it can also be seen that any competing low-redshift 
solution (shown as the dashed grey SED fits) can be ruled- 
out at high statistical significance. The fundamental reason 
for this is the additional information on the spectral slope 



long- ward of Lya provided by the UDS near-infrared pho- 
tometry. Without this additional information any sample of 
bright 2 ~ 6 LBGs candidates is vulnerable to substantial 
contamination by low-redshift interlopers, a point we will 
return to in Section 4. 



3 LYMAN ALPHA 

3.1 Lya equivalent widths 

The Lya equivalent widths were obtained by measuring the 
observed line flux, and estimating the UV continuum at Lya 
from Subaru NB921 narrow-band (NB) imaging available in 
the UDS ( |Sobral et al.|2011 |. 

The line flux was measured by subtracting any contin- 
uum from the line profile, where the continuum is estimated 
from regions either side of the line unaffected by sky line 
residuals, then integrating the flux over the pixels contribut- 
ing to the line. To ensure that the calibration of the spectra 
matches the photometry, they are convolved with the Sub- 
aru z'-band filter profile and the flux scaled to match the 
z'-band photometry (see Fig. E|. To ensure that sky resid- 
uals do not dominate the derived flux from the convolved 
spectrum, the continuum level is estimated just red-wards 
and blue-wards of Lya and is modelled by a simple step 
function. Over the small wavelength range within the filter, 
differences in true slope red-wards of Lya produce minimal 
differences to the final convolved flux. 

Absolute calibration of the spectra in this way is not ef- 
fective if the signal-to-noise is so low in the continuum that 
it is effectively measured as zero. In these cases (noted in 
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Figure 4. An illustration of the method adopted to calculate the 
Lya EWs. The spectrum is of object UUDSJ-021735. 34-051032.6 
scaled to the height of the Lya line. Over-plotted as the red solid 
line is the throughput of the Subaru z'-band filter. The spectrum 
is convolved with this filter and scaled to match the z'-band pho- 
tometry. To ensure that sky residuals do not dominate the con- 
volution the continuum is modelled as a simple step function, 
as shown by the black long dashed line, where the continuum is 
derived from regions either side of the emission lines (shown by 
the points with associated error bars). The green dot-dashed line 
shows the Subaru narrow-band NB921 filter profile. Photometry 
from this filter is used to determine the rest-frame UV flux red- 
ward of Lya emission. This is used with the calibrated line flux 
to determine the Lya EW. 



Table 1) the median correction factor measured from all the 
other spectra is used and the standard deviation of correc- 
tion factors is folded into the final error for these line fluxes. 

The final error in the Lya flux measurement includes 
pixel error estimates from the noise array, estimated error 
in the subtracted continuum measurement, as well as the 
estimated error in the derived calibration. This calibration 
error includes the error in the convolved flux, assuming er- 
rors per pixel derived from the standard deviation of pixel 
values for those contributing to continuum estimates, as well 
as the errors in the a'-band photometry. 

The NB921 photometry samples the continuum red- 
wards of Lya for z < 6.5 (see Fig. [4| and to estimate the 
continuum at Lya a flat UV continuum in F v (Fx oc A~ 2 ) 
is assumed. The NB photometry is used to estimate the UV 
continuum rather than direct measurements made from the 
spectra because sky subtraction hinders accurate continuum 
measurements in the spectra, as well as the signal-to-noise 
in any detected continuum being very low (< 2) in many 
of the spectra. The errors in the derived EWs include the 
errors in NB photometry]^] No attempt was made to correct 
either the UV continuum or Lya line flux for dust extinction 
during the EW calculation. 



3 For the two highest redshift LBGs, the Lya emission line lies 
within the extreme blue-end of the NB921 filter. Although the 
throughput of the filter is extremely low at these wavelengths, it is 
possible that the resulting EW measurements are underestimated 
by < 5%. 



Table 2. Star- formation rates derived from Lya emission-line 
luminosity and UV continuum luminosity (without any correction 
for dust extinction). The UV-derived SFR (SFRjyy) is estimated 
from the UV luminosity at A res t = 1500A using equation 1 and 
the Lya-derived SFR (SFRr^a) is estimated from the luminosity 
of the Lya emission using equation 2. 



ID 



SFR 



Lya 



SFRr 





(/M e 


yr- 1 ) 


(/M e 


yr- 1 ) 


UUDS_J021800.90-051137.8 


13.8 


± 


1.3 


30.2 


± 


2.6 


UUDSJ021616.53-050217.7 


7.4 


± 


1.3 


18.4 


± 


1.7 


UUDS.J021807.14-045841.5 


12.3 


± 


1.1 


40.6 


± 


3.8 


UUDSJ021735. 33-051032.6 


13.4 


± 


1.1 


26.3 


± 


3.6 


UUDS_J021838.90-050944.0 


16.0 


± 


1.5 


37.0 


± 


3.8 


UUDSJ021653.00-044623.3 


4.2 


± 


0.5 


30.5 


± 


2.4 


UUDS_J021922.01-045536.3 


7.5 


± 


1.0 


13.2 


± 


2.4 


UUDSJ021841. 02-051247.4 


8.8 


± 


1.5 


29.1 


± 


1.4 


UUDS.J021701. 44-050309.4 


5.0 


± 


0.9 


43.3 


± 


4.0 



Both the a'-band and narrow-band fluxes were cor- 
rected to total using measured aperture corrections from 
stars within the individual SXDS pointings. The measured 
equivalent widths were compared to estimates taken directly 
from the spectra themselves and are found to be in good 
agreement, indicating that the quoted equivalent widths are 
secure. 



3.2 Lya and UV-derived SFRs 

Both the UV continuum and Lya emission-line luminosity 
are commonly used to to derive SFR estimates. Both of these 
indicators are sensitive to the presence of dust and Lya is 
also sensitive to neutral hydrogen within the galaxy and 
the surrounding IGM. In Fig. [5] we compare the SFR esti- 
mates derived from the UV continuum luminosity using the 
|Madau et al.|(|1998| formula (corrected to a Chabrier IMF), 



with those derived from the measured Lya luminosity (e.g 
Nilsson et al.|2009 |, with no correction for dust extinction: 

L uv (/ergs -1 Hz -1 ) = 4.8 x 10 27 x SFR UV (/Meyr -1 ) 



L Lya (/ergs" 1 ) = 9.7 x 10 41 x SFR Lya (/Meyr" 1 ) 



(1) 

(2) 



The constant in equation 2 is based on the continuum lu- 
minosity at Aresi = 1500A and has been converted to a 
Chabrier IMF using a ratio of 1.65, the asymptotic ratio 
between the number of ionizing photons produced by a con- 
stant SFR with a Salpeter compared to a Chabrier IMF 



predicted by Bruzual & Chariot (2003) 



Although there is a large scatter between the two SFR 
measurements, it can immediately be seen from Fig. 5 that 
the SFR estimate provided by the Lya luminosity is system- 
atically lower than that provided by the UV continuum lumi- 
nosity, with a mean value of SFRr^a/SFRyv = 0.36 ±0.05, 
and a median of 0.4. 

We note that this is very close to the value of 0.4 which 
would be predicted by the recent study of the redshift evo- 
lution of the Lya escape f ractio n by Hayes et al. (20111. At 
low redshift Hayes et al. (20111 calculate f e Jt a by compar- 



ing the observed Lya luminosity function with a prediction 
based on the extinction corrected Ha luminosity function 
(assuming case B recombination). The lack of observed Ha 
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Figure 5. Lya-derived SFRs vs. UV-derived SFRs. The short 
dashed red line shows the mean Lya:UV SFR ratio (0.36) and 
the black long dashed line shows the ratio of 0.4 predicted from 
Hayes et al. (2011) (for f^ c a = 0.25 at z ~ 6). Objects with 
z > 6.2 are plotted as red triangles. 



luminosity functions above z ~ 2.3, means that at higher 
redshifts, the intrinsic Lya luminosity function is predicted 
from the UV, with extinction estimates being made from 
SED fitting. This calculation is essentially taking the ra- 
tio of the Lya and UV star formation rate densities, p+ 
(equations 3-5), where the star- formation rate densities are 
derived from integrating over the Lya and UV luminosity 
functions respectively, i.e. 



rLya 
J esc 



Ci 



■Obs 
P+Lya 

■Obs 
P+UV 



1+C 2 



C\ — Clv 



C 2 



-Ci 



^Lyot 



fee 



(3) 



(4) 
(5) 



Where k uv = 10.3 (Calzetti et al. 2000) and k Lya = 13.8 
and Chya = 0.445 are derived from a fit to the observed anti- 
correlation between Eb-v and f esc . Based on this formalism, 



Hayes et al. (20111 find that the escape fraction varies as 

0.25 at 



Lya 



ttr oc (1 + zj^ith £ = 2.57±g;g, giving /, 
z ~ 6, hence providing an estimate of SFRi / i /a /SFRjyv r ~ 
0.4, which agrees with the mean value derived here, within 
the errors. Consequently, based on the results of |Hayes et al.| 
(20111, we conclude that our sample of L ^ 2L* LBGs at 
z ~ 6 are consistent with a Lya escape fraction of f^c" — 
25%. 

Finally, we note that in Fig. 5, while most of the sample 



agree well with the prediction of Hayes et al. ( 2011 1, the two 



largest outliers are drawn from the high-redshift end of the 
sample at z > 6.2 and suggest a lower value of f es y c a . 



4 THE LYMAN-ALPHA EMITTER FRACTION 

In this section we use our sample of spectroscopically con- 
firmed LBGs to provide entirely new information on the Lya 



't 



* . />, -', 
. ..' 



Figure 6. A greyscale representation of the region with over- 
lapping Subaru (optical) and UKIRT (near-infrared) data in the 
UDS. Over-plotted in blue are the locations of the 20 FORS2 
masks from the UDSz programme and the area imaged by 
HST/WFC3 as part of the CANDELS survey (Grogin et al. 2011) 
is shown in red. The small black squares show all the objects in 
the field that are plausible z ^ 6 candidates following the criteria 
described in the text. The green squares show the spectroscop- 
ically confirmed z ^ 6 LAEs, while the green circles show the 
objects for which the FORS2 spectra do not provide an unam- 
biguous rcdshift. 



emitter fraction of L ^ 2L* LBGs at z ^ 6. As previously 
described, the original sample of spectroscopic candidates 
was selected using the redshift probability density function 
derived using a photometric redshift analysis (McLure et al. 
2009). However, although this strategy makes optimal use 
of the available data, in order to perform a comparison be- 
tween our results and those in the literature it is necessary 
to re-engineer our candidate selection in terms of traditional 
colour-cut criteria. Consequently, all of the parent sample 
considered for spectroscopic follow-up also satisfy the fol- 
lowing criteria: 

z < 26 (6) 

i' - z > 2.0 (7) 

z - J < 0.8 (8) 

where all magnitudes are measured within a 2" diameter 
aperture and we also require each object to be a non- 
detection in the BVR bands at the 2a level. The choice 
of a stringent i' — z ^ 2 colour cut is motivated by the fact 
that any galaxy at z 6 is predicted to display a drop of at 
least two magnitudes between the Subaru i' and z' filters as- 
suming the Madau (1995) prescription for IGM absorption. 
The z' — J ^ 0.8 colour-cut is motivated by the desire to 
effectively exclude low-redshift interlopers at z ~ 1.5, while 
maintaining sensitivity to genuine z ^ 6 galaxies with mod- 
erate reddening. 

Within the overlap region between the optical and near- 
infrared imaging in the UDS the surface density of objects 
obeying the above selection criteria is 0.023 ± 0.003 per 
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Figure 7. The evolution with redshift of the fraction of lumi- 
nous LBGs (-21.75 < M uv < -20.25) which display Lya line 
emission with rest-frame EW^ 25A. The black squares show the 
measured fractions from Stark et al. (2010, 2011) and the black 
diamond is the measured fraction at z ~ 7 estimated by Ono 
et al. (2011) from the combination of a number of spectroscopic 
samples. The red triangle shows the expected value at z ~ 7 ex- 
trapolated from the trends at lower redshift (Stark et al. 2011), 
with the fitted linear relationship and uncertainties plotted as the 
red solid and dashed lines, respectively. The blue star shows the 
results of this work based on the sample selected from the UV- 
continuum flux and with upper and lower limits as described in 
the text. 



square arcmin Therefore, within the total area covered by 
the FORS2 pointings (20 pointings, each covering an area 
of 6.8' x 6.8'), we would expect ~ 15 — 20 candidates. From 
Fig. [6] which shows the region of optical/near-infrared over- 
lap in the UDS with the positions of the FORS2 masks and 
the candidates satisfying the above criteria over-plotted, it 
can be seen that there are 19 candidates in this area, with 
14 of those having been targeted for spectroscopy. Conse- 
quently, we conclude that the sample of candidates targeted 
for spectroscopy appears to be consistent with being drawn 
randomly from the parent population of candidates satis- 
fying the colour-cut criteria described above (although see 
section 4.2.4 for further discussion). 

From our sample of fourteen objects we find seven which 
display strong Lya emission (EW^25A), which suggests a 
minimum Lya emitter fraction of 50%. However, if we ex- 
clude the AGN from our sample this fraction increases to 
54%. Either way, these are very high Lya fractions com- 
pared to previous literature results for luminous LBGs and 
clearly require some discussion. 



4.1 Comparison with the literature 

We now look to other studies of spectroscopically targeted 
LBGs at z ~ 6 to see whether this Lya fraction is com- 
parable. Stark et al. (2011) present a sample of i'-band 

> l.jj^] colour cut, 



dropouts selected with a 1775 — zsso 
which is a considerably bluer cut than that adopted here 



4 1775 and 2 8 5o are the F775W and F850LP HST/ACS filters 
respectively. 



This sample, when compared to results at lower redshift 
j Stark et al.||2010 l, provides strong evidence that the frac- 
tion of objects showing Lya in emission increases from 
2 = 4 to 2 = 6 for galaxies in the absolute magnitude 
range —20.25 < Muv < —18.75. However, in contrast, lit- 
tle evolution in the Lya fraction is seen for brighter ob- 
jects with UV luminosities more comparable to our sample 
(-21.75 < Muv < -20.25). In this luminosity range, Stark 
et al. (2011) find that the fraction of objects showing Lya 
in emission is 20% ± 8.1% and 7.4% ± 5.0% for objects with 
EWs greater than 25A and 55A respectively. 

Considering the small size of our sample, at the higher 
EW^ 55A threshold our results are in reasonably good 



agreement with Stark et al. (20111, given that one of our 
objects has EW ~ 55A (i.e. 1/9). However, at the lower 
EWJS 25A threshold our sample shows a substantially larger 
fraction of LBGs with Lya emission (~ 54%). 

4.2 Potential biases 

Understanding or reconciling the difference in the measure- 
ments made by these two studies is difficult due to the dif- 
ferent selection mechanisms used and the small numbers of 
bright objects available for study in fields smaller than the 
UDS. In fact, it is not clear that a direct comparison be- 



tween our results and those of Stark et al. (2011 \ is actually 



meaningful, given that the vast majority of our sample are 
confined to the bright half of the luminous absolute magni- 
tude bin (-21.75 < Muv < -20.25) defined by|Stark et~aT 



(20111. However, given that the trend reported by Stark et 



al. (2011), is for a decreasing fraction of strong Lya emitters 
at 2 ~ 6 with increasing UV luminosity, our results indicat- 
ing a high fraction of Lya emitters amongst luminous 2^6 
LBGs clearly requires some explanation. In this section we 
investigate several effects which could have potentially bi- 
ased our determination of the Lya emitter fraction. 

4-2.1 Reddening 

First, we consider whether our colour selection criteria are 
selecting against galaxies with internal reddening by dust 
that may produce lower Lya EWs. The i' — z' 2.0 colour 
cut does not select against reddened galaxies as this colour 
is dependent on the fraction of flux attenuated in clumpy HI 
regions in the intervening IGM, which is independent of the 
galaxy properties. Reddening tends to increase this colour 
since there is more flux at redder wavelengths still contribut- 
ing to the 2'— band flux. It is possible that the true IGM 
absorption for individual galaxies may scatter the colours of 
genuine 2^6 galaxies below the colour cut, but this scat- 
ter should be independent of the intrinsic galaxy properties. 
The z' — J colour cut may well exclude a few redder 2^6 
galaxies from our sample. This colour cut is designed to ex- 
clude low-redshift galaxies and cool galactic stars which are 
able to satisfy the i' — 2' colour cut. However, because we 
are insisting on a strict i' — 2' ^ 2.0 colour, only five addi- 
tional objects are excluded on the basis of the 2' — J ^ 0.8 
criterion over the entire UDS field, so we do not expect our 
sample to be greatly biased towards dust-free galaxies. 

However, it is important to note that if a bluer i' — 2' 
colour cut is adopted, the potential for contamination in- 
creases due to a population of massive red galaxies at 
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z ~ 1.5 with a surface density which increases significantly 
at z ^ 26. As a result, any samples selected without near- 
infrared data and a bluer i' — z colour cut will be signifi- 
cantly contaminated by low-redshift interlopers at z' ^ 26. 
Of course, until the present study, selecting significant sam- 
ples of z ~ 6 objects with z' 26 has not been possible due 
to insufficient area. 



4.2.2 Contamination from low-redshift interlopers 

To further explore how low-redshift interlopers could con- 
taminate samples of bright 2 ~ 6 LBGs, we have investi- 
gated different selection criteria using the deep HST data 
available in GOODS-S. All of our objects were found to sat- 
isfy the HST colour selection of 1775 — zg 5 o > 1.7. This lower 
colour cut is due to the difference in filter profiles between 
the HST 0850— band compared to the Subaru 2' band, as well 
as slightly increased overlap between the HST 1775 and 2350 
band filters. With the new HST/WFC3 near-infrared data 



available as part of the CANDELS survey (Grogin et al 



2011 
—j— 



Koekemoer et al. 20111 we can investigate how our 



J cut may affect these samples. We find that only two 
objects in the field satisfy 1775 — zgso > 1.7 with Z850 < 26, 
both of which are excluded by a zg 5 o — J125 < 0.75 crite- 
rion (consistent with our z 1 — J < 0.8 colour cut). However, 
in a sample selected using 2775 — 2850 > 1.3; zsso < 27, we 
find that > 90% of objects survive the zsso — J125 < 0.75 
colour cut. It is therefore clear that while a colour-cut of 
1775 — ^850 > 1.3 should be sufficient to select a clean sam- 
ple of z ~ 6 LBGs at zsso — 27, at magnitudes brighter 
than z' ^ 26 suitably deep near-infrared photometry is es- 
sential to avoid significant contamination from low-redshift 
interlopers. In fact, this point is very well illustrated by the 
SED template fits to our spectroscopically confirmed LBGs 
shown in Fig. 2. 



4-2.3 EW measurements and cosmic variance 

It is worth noting at this point that two of our objects with 
strong Lya emission have EWs that are quite close to the 
EW> 25A cut. Although our method of EW calculation is 
quite robust, not all fields have narrow band imaging just 
red- wards of Lya to supply a continuum estimate, and it is 
difficult to get high signal to noise in the continuum from 
the actual spectra at these redshifts. It is certainly feasible 
that with a different continuum estimator, or a lower limit 
derived from the noise in the spectrum, that these EWs 
could be measured to be lower, pushing down the measured 
Lya fraction. 

Moreover, it is also important to consider the problem 
of cosmic variance. Even selecting over an area of 0.25 sq. 
degrees, our sample is clearly limited by small number statis- 
tics. As illustrated by Fig. 6, previous samples of luminous 
z ~ 6 LBGs selected over the GOODS fields (the CANDELS 
imaging in the UDS covers the area of one GOODS field) 
will include very few, if any, objects as luminous as those in 
our sample. Again, this simply highlights the point that our 
sample is exploring a new area of parameter space. 



4-2.4 Photometric redshift selection 

One potential issue which needs to be addressed is the fact 
that our original object selection was performed using the 
Subaru z'-band imaging in the UDS which, for objects at 
z 6, is affected by both IGM absorption and a varying Lya 
emission line contribution. As a result, we can probably re- 
gard our estimate of 54% for the fraction of luminous 2 ~ 6 
LBGs which display Lya emission with EW^ 25A as an up- 
per limit. The fundamental reason for this is that, due to 
our original photometric redshift selection our spectroscopic 
sample is not an entirely random sampling of the parent pop- 
ulation satisfying the colour-cut criteria listed previously. In 
fact, because we targeted objects with the highest probabil- 
ity of being at 2 ^ 6, the location of the Lya emission line 
within the selection filter introduces a bias towards preferen- 
tially targeting objects with strong Lya emission. In short, 
the presence of strong Lya emission leads to an exaggerated 
Lyman break and an apparently bluer UV spectral slope, 
both of which can conspire to produce a more robust pho- 
tometric redshift solution at 2 ^ 6 than would otherwise be 
the case. 



4.3 An unbiased estimate of the Lyman-alpha 
fraction using UV continuum selection 

Given the complications introduced by selecting our original 
2^6 sample in the Subaru 2'— band, we can take advantage 
of the availability of the NB921 imaging data to investigate 
the Lya fraction in a UV-continuum selected sub-sample. In 
general, for 2 < 6.4 where the narrow-band imaging is clear 
of the Lyman-break, we expect the narrow-band fluxes to 
be the same or brighter than the 2'— band fluxes for any ob- 
ject with no Ly-a emission (this is reversed if there is strong 
Lya emission contributing to the 2'-band). A comparison 
of the 2'— band and NB921 photometry for our sample sug- 
gests that our original selection limit of 2' ^ 26 provides a 
complete sample of objects with 2921 ^ 25.7 (all objects sat- 
isfying equations 7 & 8 with 2921 ?S 25.7 also obey 2' ^ 26). 
This narrow-band magnitude cut corresponds to an absolute 
UV magnitude of 1VW = -21.0 at 2 = 6 and M uv = -21.1 
at 2 = 6.5 (assuming a flat UV slope in /„). At this limit we 
find 45% of objects (5/11) with Lya EW > 25A, confirming 
a high fraction of Lya emitters at Muv ^ —21.1. 



4.4 A lower limit to the Lyman-alpha fraction 

Using our spectroscopic sample it is also possible to estimate 
a hard lower limit to the Lya fraction amongst luminous 
2^6 LBGs, if we consider the unlikely scenario in which our 
EWJS 25A Lya detections are the only strong Lya emitters 
out of the 19 objects which could have been spectroscopi- 
cally targeted. Even in this extreme scenario the fraction of 
galaxies showing Lya with EW>25A is 37%±14% (7/19), 
where the quoted error is the poisson uncertainty. 

However, in reality, this situation is better described by 
a binomial distribution which predicts that the probability 
of finding 7 or more galaxies with EW ^ 25A from a sample 
of 14 is only p = 0.01, if the chance of "success" in each trial 
(i.e. finding EW > 25A) is p = 0.20 (as is the case for the 
most luminous 2 ~ 6 objects in Stark et al. 2011). However, 
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given the potential for bias, perhaps a more reasonable cal- 
culation is to determine the probability of finding 7 or more 
galaxies with EW ^ 25A from all 19 potential spectroscopic 
targets. This returns a probability of p = 0.07, meaning that 
we are only able to exclude the hypothesis that the true frac- 
tion of EW^ 25A Lya emitters within our sample is 20% 
at the ~ 93% confidence level. 



4.5 Dust free luminous LBGs 

Our Lya fraction results are plotted in Fig. 7 along with the 
results of Stark et al. (2010, 2011) and Ono et al. (2011). 
It is clear from this plot that our best estimate of ~ 45% 
for the Lya fraction amongst luminous LBGs at z ~ 6 is 
significantly higher than found by Stark et al. (2011) at the 
same redshift and by Ono et al. (2011) at z ~ 7 (although we 
note that, before combining their results with other studies 
taken from the literature, they find a fraction of ~ 33%). 
As far as we can tell, this result is not significantly biased 
by our sample selection and, as can be seen from Fig. 7, is 
in good agreement with a simple extrapolation of the Lya 
fraction derived by Stark et al. (2010) for luminous LBGs 
at z — 4 and 2 = 5. 

One possibility is that by sampling the bright end of 
the —21.75 < Muv < —20.25 bin, we are starting to see 
the effects of targeting the extreme tail of the population 
of the less-dusty, highly star forming systems. As we move 
to brighter and brighter objects, the extremely steep slope 
of the luminosity function would seem to require that sam- 
ples eventually become dominated by objects with low dust 
reddening and, presumably, a correspondingly high escape 
fraction of Lya photons. Although this scenario would run 
counter to the observation that, at a given redshift, brighter 
UV-selected galaxies are generally redder than their low- 
luminosity counterparts, we note that the median UV slope 
of our sample is (B) = —2.51. This is bluer than the me- 
dian value of 8 = -2m±°A 2 Q for z ~ 6, L > 0.75L* de- 



rived by Finkelstein et al. (2011 I and the variance- weighted 



mean value derived at z ~ 6.5 by McLure et al. (2011) of 
(B) — —2.05 ± 0.09, although in the absence of deeper near- 
IR photometry it is difficult to draw any firm conclusions 
regarding the UV slopes of our sample. 



4.6 A UV continuum selected sample 

As discussed previously, selecting z ^ 6 galaxy samples 
based on 2' — band photometry involves the added complica- 
tions of varying amounts of IGM absorption and Lya con- 
tamination. In this study, we were able to correct any possi- 
ble biases by using narrow-band imaging data red-wards of 
Lya to define a sample complete to a given depth in the UV 
continuum. However, to achieve this at higher redshifts it 
is clearly necessary to perform the primary sample selection 
in the near-infrared. Within this context, the new Cosmic 
Assembly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS; Koekemoer et al. 2011; Grogin et al. 2011) will 
supply deep, near-infrared, imaging (5a depth of -H160 — 27) 
over a total area of ~ 700 square arcmin. The combination 
of depth and area provided by CANDELS will allow the Lya 
fraction amongst large, unbiased, samples of z > 7 LBGs to 
be studied in the near future. 



5 CONCLUSIONS 

Targeting a sample of z ^ 6 galaxy candidates originally 
selected using a photometric redshift analysis, we obtained 
a very high redshift completeness, with 11/14 objects pro- 
viding robust spectroscopic redshifts from the detection of a 
Lya emission line. Comparing the star-formation rate esti- 
mates based on UV continuum and Lya luminosity we find 
that our sample is consistent with a Lya escape frac tion of 



fe a v c a ~ 25%, in agreement with the recent study of Hayes 



et al. (20111 



Based on our sample of L ^ 2L* LBGs at z ~ 6 we 
derive estimates for the maximum (54%±20%) and mini- 
mum (37%±14%) fraction of Lya emitters with EW> 25A. 
Testing whether these fractions are biased by our z'-band 
selection criteria, we use the available NB photometry to 
calculate the fraction of EW^ 25A Lya emitters amongst a 
complete, UV continuum selected, sub-sample. This calcu- 
lation returns a fraction of 45%±15%, consistent with our 
previous estimates. 

Our estimate of the EWJS 25A Lya emitter fraction in 
the magnitude range —21.75 < Muv < —21.25 is a factor 
of ~ 2 larger than previous estimates (i.e. 20 ± 8%; Stark 
et al. 2011). Indeed, we calculate that our sample is incon- 
sistent with a Lya emitter fraction of 20% at the ~ 93% 
confidence level. However, we also note that our sample ex- 
plores a higher redshift and luminosity range than previous 
studies. 

In summary, our results suggest that, as the epoch of 
reionization is approached, it is plausible that the Lya emit- 
ter fraction amongst luminous (L ^ 2L*) LBGs shows a 
similarly sharp increase to that observed in their lower- 
luminosity (L ^ L*) counterparts. 
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